Two-dimensional Bose gases near resonance: universal three-body effects 
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We report in this Letter the results of our investigation of 2D Bose gases beyond the dilute limit 
emphasizing the role played by three-body scattering events. We demonstrate that a competition 
between three-body attractive interactions and two-body repulsive forces results in the chemical 
potential of 2D Bose gases to exhibit a maximum at a critical scattering length beyond which these 
quantum gases possess a negative compressibility. For larger scattering lengths, the increasingly 
prominent role played by three-body attractive interactions leads to an onset instability at a second 
critical value. The three-body effects studied here are universal, fully characterized by the effective 
2D scattering length a2D (or the size of the 2D bound states) and are, in comparison to the 3D 
case, independent of three-body ultraviolet physics. We find, within our approach, the ratios of the 
contribution to the chemical potential due to three-body interactions to the one due to two-body to 
be 0.27 near the maximum of the chemical potential and 0.73 in the vicinity of the onset instability. 
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Two-dimensional quantum many-body systems have 
been, for many years, a subject of fascination for 
condensed matter and nuclear physicists alike. More 
recently, this topic also caught the attention of the 
cold atom community with the realization of quantum 
Bose gases confined to two-dimensional geometries [3- 
U These experimental studies have so far explored 
these cold atom systems at temperatures close to the 
Berezinskii-Kosterlitz-Thouless phase transition 0-0|- 
They highlighted the loss of long-range order due to 
the proliferation of vortices above the transition tem- 
perature, and the existence of two-dimensional quasi- 
condensates with algebraic long-range order and long 
wavelength thermal fluctuations below the transition. 
However, the fundamental properties of 2D Bose gases 
near absolute zero, where quantum effects are dominant, 
have yet to be addressed. In particular, on both theo- 
retical and experimental sides, very little work has been 
carried out to study 2D Bose gases near resonance. The 
main purpose of this Letter is to provide new light on the 
properties of 2D Bose gases in this limit. 

Compared to 3D Bose gases near resonance, which re- 
ceived more attention in recent years j8-ll|, 2D gases 
possess important advantages. First, the ratio between 
elastic and inelastic collision cross sections can be sig- 
nificantly enhanced when atoms are confined to two- 
dimensional traps [121 ]. Second, in 2D, trimers and few- 
body structures are all universal as the absolute energy 
scale of the spectrum is uniquely set by the two-body 
binding energy and is independent of the short distance 
property of three-body interactions This is dis- 

tinctly different from the physics of Efimov states in 3D 
as, in this case, the absolute energy scale is set by the 
ultraviolet physics of three-boson scatterings [17 1 . 

These advantages are related to the dramatic suppres- 
sion of the low energy effective interactions and phase 
shifts by coherent interference in 2D Bose gases. In fact, 
for an arbitrary repulsive interaction, the low energy two- 



body scattering phase shifts are logarithmically small in- 
dicating an asymptotically free limit. This aspect of scat- 
tering theory plays a critical role in the physics of 2D di- 
lute Bose gases. Most previous works on 2D Bose gases 
considered systems where the range of the repulsive in- 
teractions or the core size of the hardcore bosons, an, 



were much smaller than the inter-particle distances [18- 
[20j . Consequently, the results of these studies are only 
applicable when ln ^ a ^ (n is the density of bosons) is 
much smaller than unity, a limit corresponding to di- 
lute gases in 2D. Here, we focus on the physics beyond 
the dilute limit to study 2D Bose gases prepared on the 
upper branch and interacting via a resonating contact 
interaction. Such a setup can be achieved experimen- 
tally through a combination of Feshbach resonance and 
optical confinement 2ll423||. Theoretically, to study 2D 
near-resonance Bose gases, we introduce a 2D effective 
scattering length a,2D ■ This new tuning parameter is for- 
mally defined as the position of the node in the wave 
function for two scattering particles and is also identified 
as the size of the two-body bound state. In general, a^D 
can be tuned to values larger than the averaged inter- 
atomic distance and can even be infinite. 

Our study of 2D Bose gases at large scattering lengths 
unveils that near resonance the properties of these 
gases are primarily dictated by the competition between 
three-body attractive interactions and two-body repul- 
sive forces. We also show that the energetics of 2D Bose 
gases near resonance are universal as they only depend 
on the parameter na\ D . Finally, we investigate the be- 
havior of the chemical potential for a wide range of scat- 
tering lengths. We find that the chemical potential first 
increases with a-iD but very quickly reaches a maximum 
at 



ln(j: 



■ = —0.135 beyond which the Bose gas devel- 
ops a negative compressibility. Increasing a,2D further 
brings about an onset instability at ln ^ a -i — j = —0.175. 
We identify both critical values to result from the impor- 
tant role played by three-body attractive interactions. 
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Within our approach, we can estimate the contributions 
from three-body interactions to the two-body ones to be 
around 0.27 near the maximum of chemical potential and 
0.73 in the vicinity of the onset instability. 

To carry out this study of 2D Bose gases, we em- 
ploy a method previously developed to understand the 

In 
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physics of 3D Bose gases near resonance 
this approach, the chemical potential of non-condensed 
particles, and the density of condensed atoms, no, 
are first introduced as given parameters. The Hamilto- 
nian describing such a condensate interacting with non- 
condensed atoms through a short range interaction is 

H = ^(e k - v)blW + 2U n ^ 6[& k 

k k 
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Yl & I+a &t -k+a&k'+a6-k'+a + h.c. (1) 
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Here e k = fi. 2 k 2 /2m, the sum is over non-zero momentum 
states, S is the total area, and Uo is the strength of the 
bare short range interaction. Later, we will evaluate no 
and \i self-consistently as a function of the 2D scattering 
length, a,2D, and of the total density n. 

Once the full system energy density E(riQ, u) is known, 
one can calculate /x c , the chemical potential for the con- 
densed atoms, and n — no, the density of non-condensed 
atoms. This step is achieved using the following thermo- 
dynamic relations 

dE(no,n) dE(n ,n) 
A*c= tt , n = n — ; (2) 
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As hinted above, in the ground state, one requires yU c , the 
I 



chemical potential for the condensed atoms, to be equal 
to the chemical potential ii. This equilibrium condition, 
first emphasized in Ref. [26| , yields a self-consistent equa- 
tion. The evaluation of E(n , u) for a given p, and n is 



usually carried out diagrammatically [26], |2jj. To cap- 
ture the role of three-body interactions and to compare 
it with two-body contributions, we restrict ourselves to 
the virtual processes involving only two or three excited 
atoms. Truncating the Hilbert space accordingly, we can 
then sum up all connected diagrams contributing to the 
energy density. Within this truncation scheme, only the 
irreducible two- and three-body effective interaction po- 
tentials <7 2 .3 appear in the final expression for E(no,/i). 
In order to implement the self-consistency condition and 
simplify the computation of E(jiq,h), we introduce for 
the non-condensed or virtual atoms an additional pa- 
rameter n = S — /i where £(no,/z) is the self-energy. 
Physically, 77 can be understood as an energy shift due 
to the interaction between condensed and non-condensed 
atoms. Using the same series of diagrams as in our study 
of 3D Bose gases near- resonance [24j |. but carrying out 
the calculations in two spatial dimensions, we obtain for 
the energy density 

E(n ,fj.) = in 2 . 92 (2n) + ^n^Re o 3 (3n) 



with g 2 (2r\) 



H 2 4tt 
mini 3 ' 



, ff3 (3n) = 6 52 2 (2n).g 3 *(3n) 



h 2 r 

where # 3 (3n) = — / 
m J 



Aqdq G 3 (- 3rj,q) 
2n + g 2 ln 
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<7 2 .3 stand for, respectively, the renormalized two- and 
three-body interactions in a condensate. We will dis- 
cuss this point in more details below. G 3 (—377,73) repre- 
sents the three-atom off-shell scattering amplitude (cor- 
responding to the sum of all N-loop contributions with 
AT=1,2,3,...). G 3 is the solution to the following inte- 
gral equation, where h and m were exceptionally set to 
unity to improve readability, 



G 3 ("3n, /J ) 



4qdq 
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^(3T, + p>+q2)*-(pq)* K 2 V + q 



-G; (-377,9)). 



(4) 



Note that in Eqs.[3]andH B 2 = Aexp (f^J where A is 
an energy cutoff related to the effective interaction range, 
R*, via A = ^ rT . As B 2 = 32,3 are uniquely 



determined by the parameter 



ih 2 „„ „„2 



or na 



2D- 



For repulsive interactions (or positive Uo), B 2 is larger 
than A and so a 2 D is bounded from above by the interac- 
tion range R*. When Uq is infinite (hardcore potential), 
a 2 D is equal to the core size oq. For attractive interac- 



r 



tions (or negative Uq), the case we focus on here, B 2 is 
precisely the dimer binding energy, and a 2 u is the size 
of the bound state and can well exceed R* . As a con- 



sequence, na,2 D , 



the fundamental tuning parameter for 



E(no,n), can take values larger than unity. The gas can 
hence be tuned away from the dilute limit 28[. 



Before detailing our results, we need to mention, as a 
matter of completeness, that as E(no-,/J.) explicitly de- 
pends on E(no,/i), no and //, Eqs. [5] and [3] are supple- 
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Figure 1. (Color online) The chemical potential, in units of 

^jp, as a function of na\ D . The dashed (red) line is the solu- 
tion of the self-consistent equation when only two-body inter- 
actions are included. The full (blue) line is the solution when 
both two- and three-body interactions are included. This fig- 
ure highlights that the behavior of the chemical potential is 
drastically altered by three-body physics. 



Figure 2. (Color online) Three-body interaction (defined 
in Eq(3]) as a function of the energy shift r\ = E — n; 77 is 
determined self-consistently together with /j,. Inset: full and 
two-loop behavior of for small n values (respectively, full 
(blue) and dashed (red) lines). For jjL < 1, the numerical 

integration over the momentum was done from to 50 ^ B ? m ■ 



merited by the relation 



S(n 0) /x) = fx c (no,fx) + 



d In no ' 



(5) 



an extension of Hugenholtz-Pines theorem [261 . 

3D counterparts to Eqs. [2j [3] and [5] were used in 
Ref. 24| to obtain the chemical potential of 3D Bose gases 
near resonance. These self-consistent equations provided 
highly precise estimates for the chemical potential in the 
dilute limit. Near resonance, this approach predicted a 
maximum in the chemical potential and an accompanied 
onset instability. These features were fully consistent 
with the conclusions drawn from a renormalization group 
equation approach [2^]. This first study concluded that 
in 3D the dominating contribution to the chemical po- 
tential came from irreducible two-body interactions; for 
cold atoms, the three-body contribution was negligible. 
For 2D Bose gases, the story is very different: three-body 
interactions play here a much more important role as can 
be seen on Fig. [TJ 

To analyze the contribution coming from the three- 
body effect, we first solve Eqs. [2] and |3] excluding the con- 
tribution of 33, and obtain the chemical potential solely 
due to two-body interactions (see Fig. (JT|) dashed red 
line). Here, 32 is defined as the effective two-body inter- 
action renormalized by scattering events off condensed 
atoms and includes a subset of TV-body interactions de- 
fined in the vacuum [29| . Neglecting 33 interactions, the 
self-consistent equations take the simple form 



4.7T 

ln^ 



r 2a fi 
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no 
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A ln^ 



(G) 



where p, = h = 2a an d a = n al D [30(. 

solution of Eq. [6] in the limit of small a is 



n Aith 2 
m In r- 



1 =-[ln I In a I 

In — 



In47r + C] 
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(7) 



where C = In | within this self-consistent approach. This 
solution, valid in the dilute limit, agrees well with pre- 
flil . 1|1 31 j . Another solution with a ap- 
— exists in this limit but is unstable. As a 



vious studies 
proaching — 



l 2D 



2<i p 



is increased, the dilute gas solution approaches 
this higher energy unstable solution, and at the critical 
value na\ D — 1.42 x 10 -2 the two solutions coalesce into 
one. Beyond this point, no real solution to Eq. [6] exists 
revealing the presence of an instability. The basic struc- 
ture sketched here, when three-body contributions are 
neglected, is qualitatively the same as that of 3D Bose 
gases: fi is maximum when an onset instability sets in, 
and for larger na\ D develops an imaginary part implying 
the formation of molecules. 

We now turn our attention to the contribution of 
<73 (377). 33(377) is obtained by first numerically solving 
Eq. [4] for G 3 (— 3r),p) and then by carrying out the in- 
tegral involving G 3 (—3n,q) in Eq. [3] The result of this 
procedure is shown in Fig. [5] where we plot g% (3rj) . We 
chose to plot 33(377) and not 33(377) as the former is not 
cluttered by trivial effects due to g%{2rj). We identify two 
kinds of resonant scattering processes defining the basic 
structure of 33(377). The first one is a three-body reso- 
nance between three condensed atoms with zero energy 
and a dimer plus a non-condensed atom with total en- 
ergy 3?7 — £?2 • Here 3?7 is the mean-field energy shift due 
to the exchange interaction between the non-condensed 
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atom-dimer structure and the condensate. This leads 
to the first peak (from left to right) at 3rj — B^. The 
second process is a three-body resonance between three 
condensed atoms and a trimer with either binding energy 



or B [ 3 Z> (or total energies 3?7 — B^ 1 ' or 3?y 
This process produces the second and third peaks at 



? (2) 



(1) 



^ 2) ). 



3 V = B, 

? ( 2 



(1,2) 



We find numerically that B 



(i) 



1.296B 2 



and B 3 = 16.643£ 2 . These energies are fully consistent 



with the results of two previous few-body studies [13l.ll4|. 
Unlike in 3D where a logarithmically large number of Efi- 
mov states exist, in 2D there are only two trimer states. 
Remarkably, their energies are uniquely determined by 
Bi without involving an additional three-body parame- 
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ter, a fascinating feature emphasized in Refs. 

The effect of three-body scatterings on the quantum 
gas is mainly determined by the property of 173 when rj 
is relatively small. We checked numerically that in the 
limit of very small 77, 53 can be well fitted by an attractive 



interaction of the scaling form 



1m? T) ln 2 £2 ln 2 £2 

2n ■ 3tj 



captur- 



ing the dominant two- loop contribution [32j (see Fig. [2]) . 
Including the contribution due to three-body physics in 
the evaluation of the chemical potential results in two 
main effects. First, due to the attractive tail of (73 in the 
small 77 limit, as shown in Fig. [3J the instability is shifted 
away from na\ D — 1.42 x 10~ 2 and occurs at a much 
smaller value of na\ D = 3.26 x 10 -3 . At this new in- 
stability point, the chemical potential is dramatically re- 
duced, from 9.82^1 to 0.601^ when 03 is included. In 
other words, the three-body effective interaction further 
destabilizes the quantum gas. The second and equally 
important effect is that the inclusion of three-body inter- 
actions results in the appearance of a maximum in the 
chemical potential at na\ D = 0.604 x 10 -3 before the on- 
set instability occurs. The maximum value of the chem- 
ical potential is /i mQX = 1.45- iL - 2 and the condensation 
fraction at the maximum is 91%. 

Between the maximum and instability points, the 
quantum gas exhibits a negative compressibility and can 
potentially collapse into a high density phase. Although 
the fate of the Bose gases with negative compressibilities 
and the details of the corresponding dynamics are beyond 
the scope of our investigation, we speculate that in this 
regime a quantum gas eventually evolves into the droplet 
matter discussed in Ref. [Hj]. In 3D, the instability orig- 
inated from a shift of the dimers due to scatterings off 
condensates and was a precursor of the sign change of 
the effective two-body interaction (Hf. For 2D Bose 
gases, the situation is completely different. Here, the in- 
stability is a consequence of the competition between the 
repulsive two-body interaction (positive 52) and the at- 
tractive three-body interaction (negative (73) in the low 
energy limit. For a 2D Fermi gas, the Pauli blocking ef- 
fect was recently demonstrated to lead to an instability 
at a finite scattering length [23j . 

We also plot in Fig. [3] the relative weight of the three- 
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Figure 3. (Color online) Top panel: ratio between the contri- 
butions of three-body and two-body interactions as a function 
of na\ D (full red line), chemical potential for 2D Bose gases 
(full blue line). An additional metastable solution (dashed 
blue line) also exists when 53 is included. The maximum value 
of [i is 1.45^ and occurs at na\ D = 0.604 x 10~ 3 . Bottom 
left panel: condensation fraction n ( )/n as a function of na% D . 
Bottom right panel: imaginary part of the chemical potential 
when taking into account the contribution of all three-body 
recombination processes. Note that |Im fi\ <C Re /_t for all con- 
sidered na,2ij, indicating the quasi-static nature of the Bose 
gases. Hence, three-body recombination plays very little role 
in our energetic analysis and can be safely neglected for the 
range of parameters considered. 



body to two-body contributions to the chemical poten- 
tial. As anticipated, the three-body contribution is neg- 
ligible in the dilute limit when na\ D <C 1 but quickly 
becomes important as na\ D is increased. The prominent 
role played by three-body scattering leads to a maximum 
in the chemical potential before the instability point. At 
this maximum, the ratio between the three-body and 
two-body contributions reaches 0.27. The shift of the 
instability is also caused by the attractive three-body in- 
teractions. 

In conclusion, we demonstrated that the properties of 
2D Bose gases at large scattering lengths or near res- 
onance are dictated by three-body effects. We showed 
that the contributions of trimer states are universal as 
they only depend on the effective two-body scattering 
length ci2D and not on the short distance properties of 
three-body interactions; an aspect unique to 2D Bose 
gases. Our results also suggest the existence of strong 
correlations in the three-atom channel near resonance. 
This feature remains to be probed experimentally. 
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